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Abstract: ASN-Minimax double sampling plans by variables for a normally dis- 
tributed quality characteristic with unknown a and two-sided specification limits 
are introduced. These plans base on the essentially Maximum-Likelihood (ML) es- 
timator p* and the Minimum Variance Unbiased (MVU) estimator p of the fraction 
defective p. The operation characteristic (OC) of the ASN-Minimax double sam- 
pling plans is determined by using the independent random variables p*, p\ and p±, 
P2, which relate to the first and second samples, respectively. The maximum of the 
average sample number (ASN) of these plans is shown to be considerably smaller 
than the sample size of the corresponding single sampling plans. 
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1. Introduction 

For a normally distributed characteristic X ~ N([i, a 2 ) with unknown a > 
and lower and upper specification limits L and U, respectively, the fraction 
defective function p(fi, a) is denned as 

p(/i, a) := P(X < L) + P(X > U) = $ (^A + $ (^—^] , (1) 



a J \ a 



where $ denotes the standard normal distribution function. Given a large- 
sized lot and a single sample Xi, X n , (n > 3) with 



X = 

n 



n ^-^ n — 1 ^-^ 

j=i i=i 

there are two well-known procedures to determine single sampling plans. The 
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first procedure (cf. Bruhn-Suhr/Krumbholz (1990)) is based on the essen- 
tially Maximum-Likelihood (ML) estimator 



f = p( X,S) = *l^-) + z(^). (2) 



The second procedure (cf. Bruhn-Suhr/Krumbholz (1991)) is based on the 
Minimum Variance Unbiased (MVU) estimator p (cf. Kolmogorov (1953), 
Lieberman/Resnikoff (1955)). Let 



1 1 X — L \fn , 
V .. = max{0 ,---—^, (3) 

f 1 1U-X Jn 1 
W:=™ X {0,--- — ^ (4) 

*>(*):= /.^ir ^" 2 (1-^" 2 (0<x<l), (5) 

I ( ^ ^ 1 I ' 



B(x) := f b(t)dt, (6) 
Jo 

where T denotes the gamma function and B(x) and b(x) the distribution and 

n-2 

density function of the symmetrical beta distribution with parameter — - — . 
The MVU estimator p is defined by 

p:= B(V) +B(W). (7) 

The lot is accepted within the single sampling plan 

(n*, k*), if p* < k* for the first procedure, 
(h, k), if p<k for the second procedure. 

In this paper, we introduce ASN-Minimax (AM) double sampling plans for 
both procedures. For given acceptable quality level p±, rejectable quality 
level p 2 and levels a and f3 of Type-I and Type-II error, respectively, the 
AM-double sampling plans feature minimal maximal average sample num- 
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ber (ASN) while satisfying the classical two-points-condition on the operation 
characteristic (OC). We define the double sampling plans by using in the sec- 
ond stage only the information obtained from the second sample. Thus, the 
double-sampling-plan-OC can be determined by implementing the correspond- 
ing single-sampling-plan-OC. The double sampling plans are computed in a 
similar way as the single sampling plans. 

In the next section, we give two well-known theorems regarding the single- 
sampling-plan-OC for each procedure and derive with their help the corre- 
sponding double-sampling-plan-OC. The computation of the AM-double sam- 
pling plans is elucidated in the third section. In order to determine two sample 
sizes and three critical values of the two-sided AM-double sampling plan, we 
use the corresponding one-sided AM-approximation. 

2. The double-sampling-plan-OC 

Before deriving the double-sampling-plan-OC, we should bring together two 
well-known theorems from the literature on single sampling. Let 

L( n ., k * ) ( t i,<T)=P(p*<k*) (8a) 

and 

l)M = P(P < *) (8b) 

denote the single-sampling-plan-OC for the first and second procedures, re- 
spectively, and let g r be the density function of the % 2 distribution with r 
degrees of freedom. 

Theorem 1: It holds that (cf. BRUHN-SUHR/KRUMBHOLZ (1990)): 
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with 

C= ^ l)(L ~ l !l' 2 1 /''<"■/') = L + r -Ma. i,) 

w (*-. {- 

Theorem 2: It holds that (cf. BRUHN-SUHR/KRUMBHOLZ (1991)): 
-$ (-<J £ + (1 - 2B~ 1 (A;)) v /t(n-l)) |y fl _i(t)dt + 

+2j |y $' (-^ + (2 2/ -i) v /t(n-i)j • 

V«(w-l)^-i(*)dt}dy (10) 



with 



h{U -Lf 

if x + y < 1 



A(x,y)= <J Mn-l)(l-x-y) 

oo else 



S v = Vfi^—^, 5 L = Vn~^—^ and i;(y) = B-\k- B(y)). 
a a 



Definition 1: The double sampling plan by variables 



n* Zr** Zr** 
I ]_ ]_ 2 I 

A = ( ) for the first procedure, 

n 2 ^3 

A= I ~* 2 I for the second procedure, 



n 2 k. 
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with n\,n^hi,h 2 G N; ra*, n^, ni, tt-2 > 2; fc*,^,^,^,^,^^; < &2, 
k\ < k 2 , is defined as follows: 

(i) Observe a first sample of size n* for the first procedure, hi for the second 
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procedure and compute pX, pi- 
ll p\ < k*, p± < fci, accept the lot. 
If p\ > k 2 , Pi > k 2 , reject the lot. 
If A;* < p* < fcg, ki < Pi < h, go to (ii). 

(ii) Observe a second sample of size n 2 for the first procedure, h 2 for the 
second procedure and compute p 2 ,p2- 

If P2 < fcg, p2 < ^3, accept the lot. 
If P2 > fcg, p2 > ^3, reject the lot. 

Let p es * G {p*,p} denote the used p-estimator, (n, k) G {(n*, k*),(h, k)} 

the corresponding single sampling plan and A = ( 1 1 2 ], A G {A*,A} 

\n 2 k 3 J 

the corresponding double sampling plan given in Definition 1. Considering the 
independence of the samples, the OC L x (fj,, a) is given by 

L x (n, a) = P M {pT < h) + P M (pf < h, h < pf < k 2 ) 

The ASN N x (n, a) is given by 

N x (fi, a) = m + n 2 P M (h < pf < k 2 ) (12) 

with 

P(n,a)(ki < pf < k 2 ) = L(n lM )(fi, a) - L (niM) (/j,,a). 

Depending on jj, and a, L x and N x are bands in p (Figures 1, 2). The maximum 
N m axW of N x (fx,a) is for the computation of the AM-double sampling plans 
\am irrelevant. However, we calculate N max (X) as a basis for comparison with 
the ASN maximum of the AM-one-sided approximation, which is used in the 
computation of A am- A detailed description of the computation of A am is 
given in the next section. 
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3. The computation of the AM-double sampling plans 



The two-sided single sampling plan is calculated by using the one-sided ap- 
proximation (n, k), which, without loss of generality is chosen for the case of 
an upper specification limit U. For given acceptable quality level p±, rejectable 
quality level p 2 and levels a and f3 of Type-I and Type-II error, respectively, 
the one-sided approximation (h, k) is computed by fulfilling the condition 

(i) Ffl-i, %i)(0 > 1 - « 

(ii) Fn-i, 5(p 2 )(0 < (3 

(iii) ft = min (13) 



(iv) k 



$ ( —= J for the first procedure 

1 I 

B ( — I — — I for the second procedure. 

V2 2(n-l)J 



Fh-i, 6(p) denotes the distribution function of the noncentral t-distribution with 
n — 1 degrees of freedom and the non-centrality parameter 

5(p) = Vfi^ip). (14) 

Let for the case of an upper specification limit U, 4>* — \ 1 1 2 ) denote 

\n 2 h J 

the AM-double sampling plan defined by using the independent statistics 

T, = v^^V^, 2 = 1,2. (15) 

The OC and ASN for a double sampling plan (p in this case (cf. Hilbert 
(2005)) are given by 

L 4>(p) = Fni-i, Si(p)(k) + K2-1, s 2 (p)(h) (F ni -i, Si(p)(k) - ^m-i, Si(p)(h)) (16) 
and 

N 4>(p) = N (P) = ni+n 2 (F ni _ h Sl (p)(h) ~ ^m-i, 5 1 ( P )(h)), (17) 
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respectively, where 5i(p) = ^/nl <£> l (p) for % = 1,2. 
(J)* is computed by 

(i) I^(pi) > 1 - a 

(ii) L^{p 2 )<P (18) 

(iii) N max ((f)*) = minAU^ (0), 

where Z denotes the set of the double sampling plans fulfilling (18) (i) and 

(ii). The algorithm for the computation of 4>* is similar to that developed by 

Vangjeli (2009) for computing AM-double stage tests. The parameters of the 

~ / hi k\ k 2 i 
approximated AM-double sampling plan A = ~ in the two-sided 

\ h 2 k 3 J 

case are given by hi = ni, h 2 = n 2 for both procedures and 

ki = $ I -±=) ,h = ^( 4=1 ,h = * ( 1 



for the first procedure, and 

h = B (U /k 2 = B (U , k 3 = B ri 



2 2(71! -1) J' \2 2{n 1 -l)J , 6 \2 2(n 2 - I) J ' 

for the second procedure. Identically to the computation of the single sampling 
plans, we set 

L + U 

= 



0o (p) 



2 ' 
L-U 



2*-i (f ) 

and define for given p (0 < p < 1) and <r (0 < a < a (p)), /i = fi(a,p) by 
inverting (1) with /i > /i . Setting 

l a am (0;p) : = L A AM (/i(^,p),^), (19) 

^A AM (a;p) := N XAM (fi(a,p), a), (20) 
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the AM-double sampling plan Xam is given by the condition 



(i) min L XAM (a; Pl ) >l-a 

0<a<ao[p) 



(ii) max L XAM (a;p 2 ) < /? (Figure 1) 

0<cr<ao(p) 

(iii) N max ((p*) = mmN max ((p), 

4>ez 



the ASN maximum of Xam is defined 



Nmax ( A am) = max N Xam (a; p) (Figure 2) . 



(21) 



(22) 



A^m is calculated as follows: 

1. Compute the single sampling plan (n, k) and save «*,/?* 

2. Starting from a** = a*,/?** = /?*, compute A for p 1 ,p 2 ,a**, {3**, verify 
(21)(i), (ii) and vary a**,f3** until (21) (i), (ii) are true. 



Example 1 

For L — l, U = 9, pi = 0.01, p 2 = 0.06, a = /3 = 0.1, we get 

(i) (n*, fc*) = (36, 0.02645943143) and a* = 0.082, (3* = 0.1. 



a** 


13** 


A 


N max {4>*) 


minLr(cr;pi) 


maxLr(cr;p 2 ) 

(7 






25 


0.016988 


0.034411 








0.082 


0.1 








31.31538 


0.8882933643 


0.0999999893 






19 


0.029059 














26 


0.017688 


0.034554 








0.076 


0.1 








32.16417 


0.8960633610 


0.0999999882 






19 


0.029215 














26 


0.017577 


0.035291 








0.072 


0.1 








32.75441 


0.9010124424 


0.0999999889 






20 


0.029275 











where A 



AM 



26 0.017577 0.035291 
20 0.029275 



with N max (X* AM ) = 32.75439. 
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(ii) (ra, k) = (34, 0.02262119182) and a* = 0.098, /3* = 0.094. 



a** 


j3** 


A 




min Lr(o-;pi) 


maxLr(cr;p2) 


0.098 
0.098 


0.094 
0.093 


24 0.011840 0.029264 

18 0.023800 

24 0.012148 0.029093 

19 0.023424 


30.19915 
30.34629 


0.9002805676 
0.9002364848 


0.1004854552 
0.0995258042 



where \ A m 
Example 2 



24 0.012148 0.029093 
19 0.023424 



with N max {\ AM ) = 30.34628. 



For L — l, U = 9, p 1 = 0.01, p 2 = 0.03, a = (3 = 0.1, we get 

(i) (n*, k*) = (115, 0.0178762881) and a* = 0.085, p* = 0.1. 



a** 


£** 


A 


N max {4>*) 


minLr(cr;pi) 


maxLr(cr;p 2 ) 






79 


0.013777 


0.021642 








0.085 


0.1 








101.1604 


0.8948821204 


0.0999999568 






64 


0.018624 














80 


0.013902 


0.021726 








0.083 


0.1 


102.0913 


0.8972199027 


0.0999999565 






64 


0.018464 














81 


0.014029 


0.021742 








0.080 


0.1 








103.5434 


0.9008045948 


0.0999999565 






66 


0.018537 











where \* AM = 



81 0.014029 0.021742 
66 0.018537 



with N max (X AM ) = 103.5434. 
(ii) (n, k) = (113, 0.01678745123) and a* = 0.096, /5* = 0.094. 



a** 




A 




min Lr(cr;pi) 

cr * 


max Lr(cr;p 2 ) 

a A 


0.096 
0.096 


0.094 
0.092 


78 0.012471 0.020036 
62 0.017078 
78 0.012406 0.020069 
64 0.016981 


99.13899 
100.1070 


0.9000170882 
0.9000091667 


0.1014154024 
0.0993767725 



where Aa 



m 



78 0.012406 0.020069 
64 0.016981 



with N max {\ AM ) = 100.1070. 
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Remark: Numerical investigations have shown: 

(i) There are nonessential differences between N max ((f)*) and N max (XAM)- 

(ii) Since the parameters a** and f3** vary by a constant of 0.001, we can- 
not presume that the calculated Xam is the precise AM-double sampling 
plan. Finding plans that fulfill more sharply (21) (i), (ii) and feature 
an insignificantly lower N max would require an unreasonable amount of 
time. 

(iii) For the calculated double sampling plans, it holds that 

Nmaxi^Au) < N max (X* AM ) 

(Figure 3), which is consistent with the known superiority of p over p* 
from the single sampling literature. 




Figure 3: ASN bands for X* AM and Xam defined by p\ = 0.01, P2 = 0.06 and 
a = (3 = 0.1 
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